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High-precision gyroscopes are a key component of inertial navigation systems. By considering 
matter wave gyroscopes that make use of entanglement it should be possible to gain some advantages 
in terms of sensitivity, size, and resources used over unentangled optical systems. In this paper 
we consider the details of such a quantum-enhanced atom interferometry scheme based on atoms 
trapped in a carefully-chosen rotating trap. We consider all the steps: entanglement generation, 
phase imprinting, and read-out of the signal and show that quantum enhancement should be possible 
in principle. While the improvement in performance over equivalent unentangled schemes is small, 
our feasibility study opens the door to further developments and improvements. 
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I. INTRODUCTION 

High-precision inertial navigation systems (INSs) are an important enabling technology [1]. The fact they do not 
rely on external references means that they are secure against deception or jamming, which is a critical weakness of 
global navigation satellite systems. Gyroscopes and linear accelerometers are the key components of INSs since the 
time integration of their signals allows for dead reckoning, i.e. the location and orientation of an object to be known 
at all times. There is therefore considerable interest in improving their precision and performance. 

Current gyroscopes rely on mechanical rotating flywheels, vibrating microelectomechanical systems (MEMS), or 
light in ring lasers or fibre optic arrangements. The best available devices are based on ring lasers that make use of 
the Sagnac effect whereby light is sent in opposite directions around a closed loop and the two components acquire a 
relative phase dependent of the angular velocity of the loop [2]. Such a technique can measure angular velocities with 
sensitivities of around 7.8 prad/s/VHz [3], however these rely on rings with impractically large areas. The Sagnac 
effect applies equally well to atoms but the relative sensitivity to phase of matter waves is much greater than in optical 
systems. So atomic ring gyroscopes offer the potential for improved precision as well as more practical ring sizes. This 
effect has been experimentally demonstrated [4] and, while this has not yet reached the sensitivity for light, future 
experiments should be able to surpass it [5]. 

Further improvements should be possible by making use of entangled states [6] in a quantum metrology scheme. It 
is known that entanglement can enable better measurement precision to be achieved with the same resources [7-18]. 
The trade off is that entangled states are difficult to create and manipulate. This means that any realistic scheme is 
likely to involve a sequence of small entangled states and hence the improvement is likely to be modest, though this 
may still prove to be important. Quantum metrology is also known to have other advantages such as when delicate 
samples are being measured [19, 20]. 

This paper demonstrates a proof-of-principle of a quantum-enhanced atomic gyroscope. We consider the experi¬ 
mentally accessible system of an atomic Bose-Einstein condensate trapped in a rotating two-dimensional trap and 
investigate how the entangled state can be created, the rotation imprinted, and the final signal read-out, thus pro¬ 
viding an interferometric protocol used to estimate external rotations. We are able to show that by making use of 
entanglement, an improved sensitivity to rotations should be able to be achieved in this system. While the improve¬ 
ment is modest, this study shows that quantum-enhanced atomic gyroscopes are a realistic prospect and paves the 
way to further investigations that will improve the advantage and address practical issues and constraints that will 
be important in delivering a new technology. 


II. MODEL 

The physical system we consider consists of a mesoscopic sample of N bosonic atoms of mass M in an axially 
symmetric harmonic potential, with frequency ujj_ in the xy plane and Wz in the 0 axis, interacting through hard- 
core-type elastic collisions. We take fiujz to be very large compared to the interaction energy so that the dynamics 
along the z axis are frozen, i.e. all particles occupy the lowest axial energy level, thus rendering the gas effectively 
two-dimensional at sufficiently low temperatures. The trapped gas is rotated at angular frequency H around the z 
axis with the aid of an external potential which in the rotating frame appears as an anisotropic quadratic potential 
in the xy plane. The Hamiltonian in the rotating frame is given by 
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where H(<C 1) quantihes the degree of anisotropy, which is treated perturbatively in the calculations. Here, L^i is the 
angular momentum component in the z direction of the i-th atom and the —ilLz term transforms the Hamiltonian 
to the rotating frame, where Lz = is the total angular momentum of the condensate. Finally, g is the 

dimensionless interaction coupling constant which quantifies the strength of two-particle interactions and is related 
to the 3D scattering length ahy g = V^a/Xz, where Az = y^h/MuJz- Hereafter, we use dimensionless variables with 
a;_L and as units of frequency and length; consequently, energy and angular momentum are expressed in 

units of fiu}±_ and fi, respectively. 

In the isotropic case {A = 0), the Hamiltonian is exactly diagonalisable in blocks of definite total angular momentum 
for which a suitable Fock basis is given by the eigenstates of the total angular momentum | A^o, iVi,... iVt,...), where 
iVk specifies the number of atoms in the single-particle level labeled by the index k which represents a pair of quantum 
numbers (nk,m.k), the principal quantum number and the projection of the angular momentum respectively [21]. 
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When A ^ 0, the anisotropic term connects subspaces of given total angular momentum Lz separated by two units 
of h. If the anisotropy is small, a suitable basis to study the system is given by a truncated isotropic basis comprising 
the subspaces with L < Lmax; where the particular finite maximum angular momentum cut-off Lmax is chosen to 
ensure the convergence for the energies and eigenstates of the Hamiltonian [22]. In addition, due to the absence of 
any parity-breaking terms in the Hamiltonian (i.e. ones that couple even and odd angular momenta) only subspaces 
with even angular momentum parity need to be considered. 

Each subspace of definite total angular momentum can be further classified according to the concept of Landau 
levels. In the case of independent bosons, near the centrifugal limit fl « uj±, the energy levels are grouped in highly 
degenerate levels called Landau levels, roughly separated by an energy gap of 2oj±. In this sense, any basis state is said 
to pertain to the riLL-th Landau level if till = I + + (|wk| — TOk)/2] fVk- The importance of this classification 

resides in the fact that for fast rotating and weakly interacting condensates, the atomic dynamics is restricted to the 
lowest Landau level (LLL), when the standard LLL validity criterion is satisfied, i.e. the mean interaction energy is 
much smaller than the energy separation 2huj± [21]. The LLL approximation corresponds to considering only basis 
states with particles exclusively occupying single-particle orbitals Uk = 0 and mk > 0. This greatly reduces the 
Hilbert space and allows for a more computationally tractable system. However, it has been recently shown that, 
although the LLL approximation can suffy for the purpose of studying the nucleation of vortices [36] or symmetry 
breaking [22], it is necessary to go beyond this approximation in order to accurately describe the quantum states of 
the system that are useful for metrology [23]. In fact, the LLL approximation is completely unable to capture the 
possibility of creating NOON states, whereas the inclusion of a larger basis allows for it, and reveals a much richer 
system amenable to quantum metrology. In the present work, a better approximation that goes beyond the LLL is 
considered by including more Landau levels in our calculations. 

Here, it is convenient to present the stationary energy spectrum of the system, as it provides a useful insight into 
the manipulation of the ground state considered to achieve the interferometric steps. The energy spectrum for = 12 
particles obtained with a two-Landau-level basis and maximum angular momentum cut-off Lmax = fV -I- 4 is shown 
in Fig. 1. Importantly, the inclusion of a non-vanishing anisotropic term in the Hamiltonian results in an avoided 
crossing between the ground state and the first excited state near a critical frequency for any number of atoms. The 
location of this is indicated by the arrow in Fig. I. As the condensate is adiabatically brought from being at rest to a 
rotation frequency just above the first avoided crossing flniin, where the ground state contains one vortex, the system 
passes a critical frequency He where the system undergoes turbulent symmetry breaking, heralding a quantum phase 
transition. Phase transitions are well known for showing/presenting an increase of correlations between particles, 
both in the classical context and the quantum one. Particularly, in the case of quantum phase transitions, critical 
phenomena is associated with scaling of entanglement in the vicinity of the transition point, which has been found to 
be a useful metrological resource that provides an improvement in the precision of the estimation of coupling constants 
and field strengths [30-32]. 

Exact diagonalisation results show that, for even numbers of particles, the ground state at He is a strongly correlated 
entangled state well described by a two-mode approximation 

NI2 

I^Gs(f^c)) = ^C„l7V-2n)l2n), (2) 

n—0 

where the two modes correspond to the most, and second most populated single-particle states f/'i and 1 /^ 2 j respectively. 
These two states are eigenstates of the single particle density matrix (SPDM) = ('I'Gs|akai|4'Gs) and their 
populations are equal at the critical frequency, which together account for most of the population of single-particle 
orbitals. Only the LLL single-particle states having nk = 0 and Wk = 0,1, 2 take part in the expansion of f/'i and ip 2 - 
Below and near H^, ipi is a linear combination of the LLL states with mk = 0 and mk = 2, while ip 2 is proportional 
to the state with mk = 1. At He, these two most populated states abruptly interchange their composition and remain 
the same up to a certain rotation frequency greater than He. The two modes in Eq.(2) for any ground state at 
criticality that can be obtained under variations of the different parameters considered in this work, can be shown 
to be entangled by using, for example, a Von Neumann entropy criterion for mode entanglement. More importantly, 
if the same state is expressed in the atom basis, it can also be shown that any one atom is entangled with the rest 
of them. Inter-atom entanglement is responsible for quantum enhancement in metrology protocols, even when mode 
entanglement is not present [37]. 

The exact form of the two-mode approximation in Eq.(2) significantly varies depending on the relative strength 
between the anisotropic perturbation and the interactions characterised by g/A [23]. Exploration of the system for 
large values of A requires a larger angular momentum cut-off, and thus a larger basis. In this work, we restrict the 
anisotropy strength to a fixed value of A = 0.03 and varying values of g < (6/V) so that the system can be well 
described with till = 2 Landau levels and a maximum angular momentum Lmax = N + 4. In principle, however, one 
can fix the interaction strength g and vary the value of A in order to get similar results for the exact ground state 
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FIG. 1. Energy spectrum as a function of the rotation frequency for N = 12 particles as calculated with a two Landau level 
approximation. (Left) Energy spectrum for the axisymmetric case = 0). (Right) The energy levels for A = 0.03. In each 
case the interaction strength is given by gN/6 = 1 and Lmax = + 4. An avoided crossing appears in the asymmetric case at 

the location indicated by the arrow. 


form. 

When the interaction strength g is much smaller than 6/lV, the critical frequency Lie is closer to the harmonic 
trap frequency. In this case, in the absence of any anisotropy (A = 0), the lowest-lying many-body eigenstates are 
energetically close together and the ground state at the first energy crossing is quasi-degenerate involving states 
with total angular momentum D<L<N,L^1. The addition of the anisotropic perturbation explicitly breaks 
the rotational symmetry and lifts the remaining degeneracy resulting in an entangled ground state at the critical 
frequency which is well approximated by a two-mode Holland-^Burnett state, also known as a “bat” state, as shown 
in Fig.2; this “bat” state is the output state to a 50/50 beam-splitter having a twin-Fock state as an input [11]. The 
particular shape of this ground state is a direct consequence of the strong coupling of different quasi-degenerate states 
near the critical frequency connected by the anisotropic term, and the “bat-like” structure of this ground state is 
known to be robust to particle loss in interferometric schemes [24]. 

Another important feature of this “bat-like” state is related to the concept of quantum Fisher information, which 
we briefly introduce next. The quantum Fisher information determines the lowest bound for the precision of any 
number of measurements of an undetermined parameter that has been encoded into the quantum state of a system, 
independent of the measurement scheme [25]. The lowest bound is given by the Cramer-Rao inequality 
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where n denote the experimental repetitions, and Fq is the quantum Fisher information, which for a pure state is 
simply given by 

Fq = 4 [(vI/'(</)lvI/'(/,)) - l(vI/'(</)|vF(<),))|2] , (4) 

with |'I''(0)) = i9j'k(0))/i90. Although the number of measurements n is of key importance in reaching an optimal 
parameter estimation, in this paper we use the quantum Fisher information tool simply as a guide to inform the 
design of a practical scheme that shows a precision enhancement, which does not depend on the actual quantum 
Fisher information. Our parameter of interest is the rotation rate O and so we will consider the Fisher information as 
a function of O. In the case of the “bat-like” state, the quantum Fisher information, after picking up an undetermined 
linear phase in the first mode, is peaked around the critical frequency and has a relatively broad bell-like shape 
[23], which means that the unknown phase can in principle be determined with higher precision when encoded in the 
ground state of the system at any O contained in this relatively broad rotation frequency window. This is a convenient 
property since it means that an experimentalist would have a sizeable margin of error when trying to hit the critical 
frequency to prepare the initial entangled state for the interferometer. 

On the other hand, for values of the interaction strength g « 6/N, the critical frequency is farther away from the 
harmonic trap frequency, and in the absence of any anisotropy, the first ground state energy crossing is a simple 
crossing between two many-body states having total angular momentum of L = 0 and L = N. The larger value of 
the interaction strength has the effect of energetically separating the previously quasi-degenerate states so that only 
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FIG. 2. The coefficients P„ = |(A^ — 2n, 2n|^o(r2c))|^ for the two-mode approximation at the critical frequency using two 
Landau levels for = 12 particles. (Left) For this bat-like state the interaction strength is given by gN/6 = 0.44, and the 
fidelity of the two-mode approximation is |{4'o(f2c)|4'TM)|^ = 0.80, where represents the two-mode approximation to 

the state. The calculated critical frequency is Qc ~ 0.938. (Right) For the cat-like state, the interaction strength is gN/6 — 1, 
and the fidelity of the two-mode approximation is |('I'o(fIc)|4iTM)|^ = 0.70. The critical frequency in this case is flc ~ 0.823. 
The strength of the anisotropy is A = 0.03, and the truncation of the basis is Lmax = -I- 4 for both panels. 


a simple crossing remains. In this case, when a non-zero anisotropy A < 0.03 is switched on, the perturbation term 
strongly mixes only the two mentioned states at the first ground state crossing. As a result, the degeneracy is lifted 
at the simple energy crossing, turning into an avoided crossing for which the ground state is well approximated by a 
two-mode “cat-like” state shown in Fig. 2. In contrast, the Fisher information as a function of for the ground state 
in this case has a resonance-like shape, with a much larger value centred at the critical frequency, and has a width 
which is roughly two orders of magnitude smaller than the one for the “bat” state case. 

In short, for a fixed value of the anisotropy {A « 0.03), an interaction strength of g « Q/N results in a cat ground 
state at the critical frequency Lie, which exhibits a resonance-like shape for the quantum Fisher information as a 
function of LI. On the other hand, when the interaction strength has a value around g « 3/Af, the resulting ground 
state at Oc is “bat-like” shaped, and the quantum Fisher information as a function of O has a broad bell-like shape 
centred around Oc [23]. This distinctive behavior of the system under changes in the interaction strength makes the 
rotating condensate a convenient system to create and tune a variety of different entangled ground states which can 
be used as initial states for quantum metrology schemes. In this article, we propose an interferometric scheme that 
uses the rotating condensate to measure external rotations, and we show that it can achieve sub-shot noise precision 
for all the different types of entangled states that can be produced by varying the interaction strength, which in turn 
can be achieved using Feshbach resonances [26]. The interferometric scheme is conceptually simple and within reach 
of current technologies. 


III. QUANTUM INTERFEROMETRIC SCHEME 

The interferometer is implemented as a quantum metrology protocol which can be divided into four main stages: 
adiabatic preparation of the initial entangled state, non-adiabatic coupling of the rotating condensate to the external 
system whose rotation frequency is to be measured, acquisition of an internal phase which encodes the external 
rotation into the ground state of the rotating condensate, and a read-out stage. Here, we show that after the phase 
acquisition, practical read-out schemes can determine the external rotation with a precision better than the classical 
(shot-noise) value for the parameters studied. 

The interferometric protocol requires the condensate rotation frequency to be varied with time. For simplicity, we 
assume a linear dependence with time as 

f2(t) = Qo + it, (5) 

where 7 is the constant rate of change for the rotation frequency, and Qq is a particular initial rotation frequency. 
This linear dependence is the most simple and tractable way of modelling the dynamics of the system, however, it is 
flexible enough since any particular functional dependence of r 2 (t) with time can, in principle, be approximated by 
a series of different linear passages. For the purpose of this work, we choose the simplest linear dependence in each 
stage in order to obtain a proof-of-principle result that exhibits sub-shot-noise precision, and thus, there is plenty of 
room for optimisation in this regard. In experiments, the variation of Q(t) as it appears in Eq.(5) corresponds to the 
variation of rotation frequency of the anisotropic potential ordinarly achieved by a stirring laser beam that rotates at 
the rate Q(t) in the laboratory frame, and has a fixed anisotropic profile, 2AMuj']_{x'j —yf) in our case, in the rotating 
frame [34? , 35]. 

In order to simulate the dynamics of the system, we follow [22] and solve the time-dependent Schrodinger equation 
id/dt\'i/{t)) = H(t)), where the only time dependency in H{t) comes from the term —Vl.{t)Lz. Since in the absence 
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FIG. 3. Predicted adiabatic evolution time as a function of the interaction strength for different numbers of particles. The full 
evolution from = 0.4 to = Sic is divided in smaller linear segments with different angular acceleration rates as given by 
Eq.(7), and a Hxed increase in rotation frequency of 1% of the trapping frequency. Results are shown for a trapping frequency 
of uj± = 27 r X 2.1 kHz. 


of anisotropy (A = 0), the Hamiltonian is exactly diagonalizable in blocks of definite total angular momentum with 
eigenvectors [23], we expand |lE'(t)) in the many-body angular momentum basis as |'F(t)) = and, by 

projecting the time-dependent Schrodinger equation onto the state |$j), we obtain a first order system of differential 
equations for Cj(t), 


i 

subject to the initial condition |'I>(t = 0)) = = 0)|$i), which we solve numerically using a Fehlberg fourth- 

fifth order Runge-Kutta method. Due to computational constraints, all the numerical calculations in this article are 
restricted to fV < 12 atoms, two Landau levels and a basis containing angular momentum subspaces of L < TV -|- 4. 
The details of each stage in the interferometric protocol are discussed in the following sections. 


A. Adiabatic preparation of the initial entangled state 


The first stage of the protocol is the creation of the initial entangled state at the critical frequency He. For 
this purpose, we start with a non-rotating condensate in the ground state, and adiabatically ramp up the rotation 
frequency until the condensate reaches the critical frequency. In experiments, the gas can be initially cooled down 
and condensed into the ground state at a moderate rotation rate Hg for which the energy gap to the first excited 
state is large enough to allow full condensation into the ground state by relaxation of the quantum gas [22]. Here, we 
assume this to be the case and consider nucleation of the entangled state by starting at Hg = 0.4 and adiabatically 
increasing the rotation frequency. 

The main challenge of adiabatically nucleating the entangled state is to ramp up the rotation frequency from Hg, 
slowly enough so that the system always remains in the ground state during the evolution, but fast enough so that 
the total evolution time is shorter than the average condensate lifetime. In order to estimate the minimum time 
needed to nucleate the entangled state, the full evolution from Hg to He is divided into n smaller linear segments 
with fixed rotation frequency increments i5H, and corresponding rotation frequency change rates 7 ^ which give the 
fastest evolution time for the i-th segment while preserving the adiabatic condition. For numerical calculations, we 
use 5H = 0.01 and the maximum acceleration rates 7 ^ are estimated using 


af;2. 
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where AE^i^ is the smallest energy gap between the ground state and the next excited state for the z-th segment, and 
Poi is the maximum tolerance for the transition probability between these two eigenstates after the adiabatic passage 
[22], whose value we take to be poi = 0-01 for numerical calculations. In fact, Eq.(7) is known to be quite restrictive as 
an estimate of the maximum acceleration rate consistent with the adiabatic condition [27], but it provides a relatively 
easy way to estimate adiabatic evolution times. 

We show the calculated minimum evolution time estimated using numerical results from exact diagonalisation for 
AEmin, and combining Eq.(5) and Eq.(7) in Fig.(3) as a function of the interaction strength. A number of important 
conclusions can be drawn from these results. Firstly, nucleation of the entangled state via this adiabatic process is 
clearly limited to condensates containing small numbers of atoms, due to the fact that typical condensate lifetimes are 
of the order of 16 s [28]. Thus, our scheme could be implemented in optical lattices of locally anisotropic rotating tight 
potentials (wj, « 27r x 2.1 kHz) with thousands of sites with mean density between N ^ 5 and « 10 atoms/site, 
which have been demonstrated in experiments [29]. Secondly, even for small numbers of atoms in relatively tight traps, 
the adiabatic evolution is feasible only for values of the interaction strength of g < (6/lV) x 0.6, which corresponds to 
the nucleation of bat-like entangled states. In contrast, practical adiabatic nucleation of cat-like states would require 
much tighter traps or great improvement in condensate lifetimes. Therefore, we mainly focus on the case of a bat-like 
state as the initial entangled state for our interferometer in the rest of the article. In actuality, this is not much 
of a limitation in quantum metrology schemes since it has been shown that bat-like states have the same potential 
of achieving nearly Heisenberg-limited precision as cat-like states in other geometries, even outperforming the latter 
when particle losses are considered [24]. Moreover, for the particular geometry considered here, it has also been shown 
that a bat-like setup gives a more sizable margin of error in producing the entangled state, as opposed to the cat-like 
case [23]. 

For notational purposes, we represent the operation performed during this first stage as 

|4-i(H = Oe)) = HA|T = 0), (8) 

where Ua represents the unitary operator corresponding to the adiabatic evolution, which transforms the non-rotating 
ground state |L = 0) to the entangled ground state at the critical rotation frequency. 


B. Coupling to an external rotation 

The next stage is to expose the system to an external rotation that we wish to measure. This shifts the rotation 
frequency of the condensate in a rapid non-adiabatic way and takes it to a regime where its natural decomposition 
into eigenstates of the Hamiltonian at the shifted frequency allows the acquisition of an internal phase by means of 
simple free time evolution. This internal phase encodes information about the external rotation into the quantum 
state of the condensate, and the details of this acquisition are considered in the next stage of the protocol. 

The external rotation coupling can be modelled as having the following effect on the quantum state of the system 

= f^A)) = f^extl'l'l(f^c)) = ^ Ci(HA)|4'i(HA)), (9) 


where flgxt is the external rotation we wish to measure, and Ha = He — Hgxt- Also, l'I'i(HA)) is the z-th lowest 
lying eigenstate in the stationary energy spectrum of the Hamiltonian at H = Ha, and ^(Ha) = (4'^ (Ha)14*1 (He)) 
are the expansion coefficients. We have made two important assumptions in Eq.(9). The most important one is 
the assumption of very fast non-adiabatic evolution during the coupling to the external rotation, so as to obtain a 
high fidelity |(4>n(HA)|'I'i(Hc))|^ « 1 between the initial entangled state and the final state after the coupling. This 
condition ensures that we retain the entanglement generated in the first stage before subjecting the system to a free 
time evolution. In fact, this is nothing but the main premise of enhanced quantum metrology, i.e., the encoding of 
unknown parameters that we wish to measure into an initial entangled state. In addition, this assumption provides 
a simple picture where the encoding of any information about the external rotation occurs only through the next 
stage in the interferometric protocol; that no information about the external rotation is gained during this sudden 
coupling is a consequence of 9|4'n(HA))/9Hext ~ 0, resulting in a null Fisher information of the final state after the 
rotation frequency shift. Incidentally, this assumption also allows for shorter numerical simulation times. The second 
assumption consists in considering the external rotation to be small enough compared to the critical frequency, so 
that the resulting many-body state after the coupling can be described by a two-level superposition involving only 
the ground state and the next excited state at Ha for a wide range of angular acceleration rates at which the coupling 
takes place. This is not a crucial assumption, since our numerical codes are not limited to this two-level approximation 
and our scheme does not rely on it. However, it will simplify the analysis in Section HID. 



The sudden condition for the external rotation coupling imposes a limitation on the maximum value of the external 
rotation |f2ext| that can be considered for a given maximum angular acceleration rate of the frequency shift. An 
estimate of the maximum value of the external rotation consistent with the assumption of non-adiabatic evolution 
can be obtained from the condition of the sudden approximation [27] 


T< 


1 


( 10 ) 


where (•) denotes expectation values with respect to the initial state |'I'i(r2c)), and T is the total time spent in the 
evolution, and (Ai7)^ = {H ) — (i7)^ is given in terms of the time-averaged Hamiltonian 

H=^£H{t)dt. ( 11 ) 

Combining these two expressions with Eq.(l), the sudden coupling assumption implies 


If^extl < 



( 12 ) 


where 7max is the maximum angular acceleration rate of the frequency shift, and (AL^) is the expectation value of 
the standard deviation of the total angular momentum with respect to |'I'i(r2c)). Since the initial entangled state 
(either bat-like or cat-like) consists of a superposition of angular momentum eigenstates with L < N, as our numerical 
calculations show, {ALz) is at most of order N. Therefore, a larger number of particles further restricts the maximum 
value of the external rotation that is consistent with the sudden assumption. Nevertheless, since we have establi shed 
that the hrst stage of the protocol already limits the feasibility of the scheme to small numbers of atoms, the -y/l/A 
factor in Eq.(12) represents only a variation of less than one order of magnitude for [Hextl- On the other hand, 
7max has a broader range of possible values. Here, we assume the external rotation to be less than one percent of 
the trap frequency in all calculations, for which Eq.(12) gives a maximum angular acceleration rate of the order of 
7max ~ 0.50 X 10“^. In fact, our numerical calculations show that a hdelity greater than 0.96 can be achieved between 
the initial and Hnal state for the external rotation coupling operation with TV < 12 at this given maximum angular 
acceleration rate, and jOext] < 0.24 x 10“^. 


C. Phase imprinting 

The next stage of the interferometric protocol consists of a simple free time evolution at the shifted rotation frequency 
Ha, for a time r. The effect of this operation is that of creating a relative phase between the two components |5 'o(Ha)) 
and |4 'i(Ha)) that appear in the expansion of I'I'ii(Ha)) in Eq.(9). This phase encodes information about the external 
rotation into the quantum state of the rotating condensate. 

The phase imprinting operation can be expressed as 

|^iii(Ha;t)) = 17te(Ha;t)|«'ii(Ha)) = ^ c,(HA)e-*^®‘(''^)n^*(f^A)), (13) 


where 17te(Ha; t) = exp (—ti7(HA)T) is the time evolution unitary operator, and AAi(HA) is the energy gap between 
the ground state and the i-th excited state at Ha- At this point, the value of the external rotation has been encoded 
into the quantum state of the condensate and it is ready to undergo a measurement stage. 


D. External rotation read-out 

The read-out procedure is the final step and consists of undoing all the previous steps before the application of 
the phase. This is similar to the role of the second beam splitter in a standard Mach-Zehnder interferometer. After 
this, a suitable measurement is performed on the system to determine the external rotation. The particular choice of 
measurement scheme determines the precision with which the external rotation is estimated. Finding a practical way 
of optimising the measurement scheme is a rather hard problem, and there is no general strategy for it. However, a 
good starting point is to calculate the best precision possible using the quantum Fisher information, which provides 
a benchmark to aim for. 
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The fact that the quantum state of the system in Eq.(13) depends on f^A in a rather non-trivial way, implies 
that there is no simple analytical expression for the quantum Fisher information of this state in general. However, a 
compact and insightful result can be obtained under some reasonable simplifying assumptions. Although our numerical 
simulations are not restricted to these assumptions, it is useful to have an analytical expression that provides additional 
insight and possibly serve as a guide to inform the design of the read-out stage. The first assumption consists of 
considering the waiting time r to be small, so that we can neglect any terms of order O(t^) and higher. Secondly, we 
assume that at any point of the interferometric protocol, the population of excited states higher than the ground state 
and first excited state is negligible, thus making it possible to think of the system as a two-level one. This assumption 
is equivalent to considering only small external rotations, since the coupling to such small external rotations implies 
that, after the coupling, the quantum state of the system can be accurately reproduced by the apropriate superposition 
of the two lowest eigenstates of the Hamiltonian at the shifted rotational frequency. Under these assumptions, the 
quantum Fisher information right after the phase acquisition stage can be shown to be approximated by 

Fq [|'I'i„(Ha;t))] « 4r"(vI/i(He)|(AL)2|d/i(H,)). (14) 

This is a very interesting result. It implies that, as long as the waiting time and the external rotation are small enough, 
the lowest bound for the precision of rotation measurements is completely determined by the angular momentum 
fluctuations of the initial entangled state. A numerical calculation of the fluctuations of angular momentum shows 
that for the cat state, as well as for the bat state and intermediate states, ((AL)^) oc demonstrating that at least 
in this regime, the system is capable of delivering Heisenberg-limited scaling. As the waiting time starts to increase, 
going beyond the approximation used to obtain Eq.(14), our numerical calculations suggest that the quantum Fisher 
information departs from the parabolic profile and shows oscillatory behaviour with respect to r, which is heavily 
dependent on the external rotation and particularly, on the number of atoms. Consequently, the scaling with N 
beyond small waiting times becomes rather irregular. Therefore, finding the optimal waiting time involves calculating 
the longest waiting time that provides a better signal, but does not cross over to the region of irregular scaling with 
N. An analysis of the optimal measurement is beyond the scope of the present work. 

This result for the quantum Fisher information suggests that we look for a measurement of total angular momentum 
in a basis where the hnal output of the interferometer has a large variance in angular momentum in order to get close 
to the lower bound given by Fq.(14). 

One possible way to accomplish this measurement scheme is the following. After waiting for a time r, the condensate 
is quickly decoupled from the external rotation, meaning that the condensate returns to the rotation rate He under 
the sudden approximation. At this point, the two lowest-lying energy levels have the largest population. Then, we 
adiabatically switch off the anisotropy, thus bringing the system to a basis where the eigenstates of the Hamiltonian 
are also eigenstates of the total angular momentum with the ground state having L = N and the hrst excited state 
L = N — k, where the value of k depends on the magnitude of the interaction strength. For gN/6 > 0.6 the value 
oi k = N, and for smaller values of g, we have k = 1. After the adiabatic anisotropy switch-off, we end up with a 
condensate occupying mostly the first two lowest-lying energy levels, which in the case of gN/6 > 0.6 means that this 
final output state right before an angular momentum measurement has a large variance in L. 

We now proceed to calculate the external rotation precision obtained from a total angular momentum measurement 
after bringing the system to the measurement basis mentioned above. In this case, it suffices to consider projective 
measurements corresponding to the projective measurement operators E{L) = where Idti) is the f-th energy 

eigenstate at Hj. when A = 0, having L\'ifi) = and thus L labels the results of the measurement, in this case, 

a well-defined total angular momentum. The probability density for result L, given the external rotation Hext, is 

p(L|Hext) = (^iv(f^c;T)|F(L)|4-iv(H,;r)), (15) 

where |4tiv(Hc;T)) is the final output state to the interferometer protocol right after the adiabatic anisotropy switch- 
off. 

Consider now n measurements of total angular momentum with results Li,..., which are randomly distributed 
according to the distribution in Fq.(15). The external rotation is estimated using the function Lest = Yl^=iFi/n, 
which is nothing but the mean value of the total angular momentum of the sample. This function is called the 
estimator function, and how precisely the n measurements are able to determine the external rotation depends on its 
particular functional form. Before the external rotation can be estimated using this function, it needs to be “corrected” 
to account for the difference in units [25]. Thus, the external rotation estimation given by the units-corrected mean 
value of the sample distribution of the estimator is 

(f^est) = | , (16) 

|u(Aest)/uHext| 

where (Test) is the mean value of the sample distribution of the estimator function Lgst- It is worth emphasising 
that this is not a quantum mechanical expectation value, but a statistical one. Nevertheless, due to the fact that 
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FIG. 4. Standard deviation for the distribution determining the external rotation as a function of the same rotation. Here, 
the estimator used to build the distribution is the total angular momentum of the interferometer output. In addition, the 
standard deviation has been scaled by a factor equal to the waiting time, in this case r = 10, in order to benchmark against 
the shot-noise limit for unentangled particles. Also, gN/6 = 0.44 and A — 0.03. 


the estimator functional form is an arithmetic mean, we can invoke the central limit theorem for large samples with 
n ^ 1, in this case 


(Lest) = (L), (17) 

where the right-hand side average is a quantum expectation value giving the mean value of the “population” distri¬ 
bution assumed within the central limit theorem context. In the same way, the variance of the sample distribution 
can be taken as 


((ALest)") = ((AL)2)/n, (18) 

according to the central limit theorem. Consequently, the precision with which the external rotation can be estimated 
with this measurement scheme is given by the error-propagation formula 


Allext = \/((AfI ext)^) 


yLl |d(Lest)/^^ext| 


(19) 


We show the precision obtained from a full simulation of the interferometric protocol for a particular set of param¬ 
eters in Fig.(4), where sub-shot noise precision is achieved for the measurement scheme discussed above. In Fig.(4), 
we have factored out ^/n appearing in Eq.l9, as well as which are both independent of the number of particles, 
in order to directly compare the obtained precision with the shot-noise limit l/y/N. This means that the precision 
shown in the mentioned figure above can be improved by a factor corresponding to the evolution time r and Ijy/n, if 
the protocol is repeated n times. Here, it is worth mentioning that any given estimator generally performs differently 
for different values of the unknown parameter, and this can be seen in our calculation corresponding to Fig. (4). In 
particular, the divergence at He is a direct consequence of a finite numerator and a denominator that goes to zero in 
Eq.(19). The error-propagation formula is known to result in divergences like this for particular values of the unknown 
parameter. However, these divergences are not fundamental and can be avoided in different ways, for instance, by 
implementing a deliberate bias in the unknown parameter or by using a different estimation strategy, such as Bayesian 
estimation. Although we have only presented results for a frequency window about ±0.5% of uj± around the critical 
frequency, our simulations show that sub-shot-noise scaling remains up to external rotation values in a frequency 
window ±5% of uj±_ wide. However, meeting the criterion for the sudden approximation for larger external rotations 
implies a larger rotational frequency acceleration, which might not be realistic. Additionally, we have also observed 
similar sub-shot-noise scaling for larger values of the interaction strength 1 > gN > 0.44. 
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FIG. 5. Standard deviation for the distribution determining the external rotation as a function of the same rotation. Here, the 
estimator used to build the distribution is the population of the ground state of the interferometer output. In addition, the 
standard deviation has been scaled by a factor equal to the waiting time, in this case r = 10, in order to benchmark against 
the shot-noise limit for unentangled particles. Also, gN/6 = 0.44 and A — 0.03. 


Although the measurement strategy mentioned above can deliver sub-shot noise precision, it presents a serious 
challenge for experimentalists, due to the adiabatic anisotropy switch-off. In contrast to the adiabatic nucleation of 
the entangled state, we need to preserve the population of several low-lying energy eigenstates, since the state of 
the condensate right after the decoupling from the external rotation can contain non-negligible overlaps with higher 
excited states, particularly for larger values of the waiting time r. Since very short waiting times are not suitable 
for practical reasons, it is advantageous to find an alternative measurement scheme. An obvious choice is a binomial 
estimation, where the two outcomes correspond to the probability of finding the system in the ground state or in any 
other excited state, regardless of which one exactly. In this case, the approach is exactly the same as before, but this 
time, we have only two measurement operators 

4(a; = 0) = |^o)(^o|, (20) 

and 

Ai(x = l) = l-|§o)(^o|, (21) 

where a; is a random variable representing the outcome of a measurement in this second scheme, which can be a; = 0 if 
the condensate is found in the ground state, or a; = 1 otherwise. Taking the estimator function as an arithemetic mean 
as before, Aest = again invoke the central limit theorem and calculate the units-corrected standard 

deviation of the sample distribution for flest corresponding to this second approach. Notably, the precision obtained 
with this approach is slightly better than the first one but, more importantly, is more amenable to experiments. Full 
simulation results of the interferometric protocol using this second approach are shown in Fig. (5) for the same set of 
parameters as in Fig. (4). Although we have not pursued a detailed study of the dependence of the interferometric 
protocol on the waiting time, it is worth mentioning that increasing waiting times can have a significant impact on 
the precision obtained from the protocol, as it is illustrated by the particular example of Fig. (6). 


IV. CONCLUSIONS 

In this paper, we have presented a proposal for an interferometric protocol to measure external rotations using a 
rotating anisotropic Bose-Einstein condensate. For simplicity we have assumed that particle losses are rare and can 
be neglected. Future work will consider the effect of losses, which connect the odd and even number subspaces and 
so complicate the dynamics. The protocol we have presented should be able to deliver sub-shot noise measurements 
as demonstrated by the quantum Fisher information of the system after an external rotation has been applied to 
















12 



Q 

_ i_ _ __1__ _i _ _i _1 

VT /(T /¥ ,Jlo /i2 

• N=n ■ Ar=10 * A^=8 o n=6 □ N=4 

FIG. 6. Standard deviation for the distribution determining the external rotation as a function of the same rotation. Here, the 
estimator used to build the distribution is the population of the ground state of the interferometer output. In addition, the 
standard deviation has been scaled by a factor equal to the waiting time, in this case r = 50, in order to benchmark against 
the shot-noise limit for unentangled particles. Also, gN/6 = 0.44 and A — 0.03. 


it. Importantly, the feasibility of sub-shot noise measurements using two different measurement schemes was demon¬ 
strated for a particular set-up, which is a real strength of the presented protocol. Out of these two measurement 
schemes, the one that delivers the best precision is a simple projective measurement which only requires measuring the 
population of the ground state of the system at the output to the interferometer. Although it has not been presented 
in this study, a Bayesian estimation strategy can also be used to deliver the same results presented in this paper. 
The Bayesian strategy has the advantage of allowing any prior knowledge of the unknown external rotation to be 
accounted for, such as any previous measurement of the external rotation. Additionally, the interferometer is highly 
tunable, including the form of the entangled state generated inside of it, and it is amenable to a proof-of-principle 
experiment. Finally, in this investigation, we have barely scratched the surface of what the system could be capable 
of doing, as a full exploration of all its parameters is still under current investigation. 
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